Low-energy effective field theory for 
finite-temperature relativistic superfluids 
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We derive the low-energy effective action governing the infrared dynamics of relativistic superfluids 
at finite temperature. We organize our derivation in an effective field theory fashion — purely in terms 
of infrared degrees of freedom and symmetries. Our degrees of freedom are the superfluid phase 
ip, and the comoving coordinates for the volume elements of the normal fluid component. The 
presence of two sound modes follows straightforwardly from Taylor-expanding the action at second 
order in small perturbations. We match our description to more conventional hydrodynamical ones, 
thus linking the functional form of our Lagrangian to the equation of state, which we assume as 
an input. We re-derive in our language some standard properties of relativistic superfluids in the 
high-temperature and low-temperature limits. As an illustration of the efficiency of our methods, 
we compute the cross-section for a sound wave (of either type) scattering off a superfluid vortex at 
temperatures right beneath the critical one. 



I. INTRODUCTION 

The low-energy, long-distance dynamics of superfluids 
are usually studied at the level of hydrodynamical equa- 
tions. At finite temperature, these describe two fluids 
interacting in a non-trivial way. For instance, they can 
exchange charge and mass on top of energy and momen- 
tum. The standard derivation of the equations of motion 
and of the so-called constitutive relations is somewhat 
cumbersome, requires a clever use of galilean invariance, 
and an extensive one of thermodynamics (see e.g. [H, 0]). 
The generalization to relativistic superfluids is possible 
but not particularly transparent. 

Here, we present an alternative approach based on ef- 
fective field theory. We extend previous work on solids 
and ordinary fluids [§l-fl0| to the case of (non-dissipative) 
relativistic superfluids. We see several advantages in our 
approach: 

• The system is described by a field theory, in terms 
of a local action. The action encodes all the dy- 
namics in the most compact way, as usual. The 
variational principle is the standard one for local 
field theories, with no Lagrange multipliers nor con- 
straints nor additional restrictions on the fields' 
variations. 

• The starting point is the long-distance degrees of 
freedom and the symmetries acting on them. The 
action is the most general one compatible with the 
symmetries, organized as a derivative expansion. 
In particular, we will not try to cook up an ad-hoc 
action that reproduces the known hydrodynamical 
equations for superfluids. We will instead derive 
the same equations straightforwardly and system- 



atically from a different set of principles- 
effective field theory. 
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• Implementing Poincare invariance is straightfor- 
ward: the local fields that parametrize the infrared 
degrees of freedom are Lorentz scalars. Laundau's 
clever use of galilean invariance is here replaced by 
the more prosaic "contract the indices" rule. 

• Thermodynamics is not needed. More precisely: 
Since the dynamics descend from an action, they 
are non-dissipative by construction. In such a case, 
thermodynamics is needed just to establish a dic- 
tionary between our field variables and the stan- 
dard hydrodynamical and thermodynamical ones 
(sect. IIVj) . Once this is done, one can just forget 
about the thermodynamical interpretation of the 
system and use the action to compute observables. 

Besides the physical importance of relativistic effects 
for certain superfluid systems — e.g. for neutron stars — 
keeping the analysis fully relativistic from the start 
presents technical and conceptual advantages. At the 
technical level, it is just simpler and neater than the 
non-relativistic one — all our dynamical fields transform 
linearly in an irreducible representation of the Poincare 
group (the scalar one), so that enforcing Poincare in- 
variance is trivial. By contrast, for instance, the velocity 
fields of the standard non-relativistic hydrodynamical de- 
scription transform non-linearly under Galileo boosts. At 
the conceptual level, the relativistic analysis dodges cer- 
tain somewhat misleading accidental degeneracies that 
arise in the non-relativistic limit but that are not enforced 
by the symmetries. For instance, in the non-relativistic 
limit mass density and mass current are proportional to 
charge density and charge current. When relativistic ef- 
fects are included, the degeneracy is gone: 'mass' is not 
conserved anymore — it is not even well defined — while 
charge is, as enforced by a U(l) symmetry. The non- 
relativistic limit can be taken straightforwardly when 
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(and if) needed. 

For genuinely hydrodynamical questions our approach 
is probably no better than the standard one. If, for in- 
stance, the goal is to solve the hydrodynamical equations 
with certain boundary conditions, then the fact that we 
have a neater way to get to those equations is of no use — 
we still have to solve them! There are however certain 
questions for which our formalism is manifestly more con- 
venient. For instance, perturbation theory is straightfor- 
ward and systematic, as usual for a local field theory with 
an action: In sect. Illll we show that there are two propa- 
gating sound modes. In sect. IVIII we consider how these 
interact with a superfluid vortex. Computing the associ- 
ated scattering cross-sections is straightforward, via the 
usual Feynman rules. 

From the field-theoretical viewpoint, this paper builds 
on previous work on the field theory of fluids |8l4l0| and of 
zero-temperature superfluids 111], and carries similarities 
to ref. [12j . where supersolids are discussed in effective 
field theory terms. From the hydrodynamical viewpoint, 
our analysis will mirror that of [H El, with a difierent 
starting point. 

II. THE SETUP 

We begin by briefly reviewing the field theoretical de- 
scription of a zero-temperature superfluid (see e.g. [U). 
In field-theory terms, a superfluid can be thought of as a 
system carrying a conserved 17(1) charge in a state that 
(i) has finite density for this charge and that (ii) sponta- 
neously breaks the corresponding U(l) symmetry. The 
spontaneous breaking — or Bose-Einstein condensation — 
is in a sense physically caused by having finite charge 
density, but, once we consider a state that has both, we 
can discuss their implications in any order |15| . We find it 
more convenient to start with the spontaneous breaking. 
It implies the existence of a gapless excitation ip — the 
Goldstone boson — which non-linearly realizes the U(l) 
symmetry: 

ip —> ip + a . a — const (1) 

Barring accidents, one can assume that the low-energy 
dynamics involve just this degree of freedom. The low- 
energy Lagrangian for ip should then be the most general 
one compatible with the shift-symmetry (fT]) as well with 
Poincare invariance, organized as a derivative expansion. 
At lowest order in derivatives, it takes the form 

C = P(X) , X = dpW1> > (2) 

where P is, for the moment, a generic function. The 
current associated with the U{1) symmetry (JTJ) is 

f = 2P'{X)d fl ip . (3) 

We see that for the superfluid state to have finite charge 
density, one needs finite ip. In a state of uniform charge 



density and vanishing spacial current, 

iP = fit , (4) 

where fi is the chemical potential The function 

P(X) turns out to be nothing but the equation of state, 
interpreted as a relation between the pressure p and 
the chemical potential, p = P{—fi 2 ) [Tl|. Notice that 
from the field-theoretical viewpoint, the system's be- 
having as a fluid — that is, its obeying the equations of 
hydrodynamics — is quite non-trivial, and surprising. The 
corresponding fluid four- velocity field is d^ip, suitably 
normalized. In particular, it obeys a relativistic gener- 
alization of the irrationality condition, i.e., it describes 
potential flow. 

Now, Laundau's intuition is that a superfluid at fi- 
nite temperature should be describable as a mixture of 
two mutually interacting fluids, one that behaves like the 
zero temperature superfluid we have just described, with 
irrotational flow only, and one that behaves like a nor- 
mal fluid. The intuition is substantiated by the following 
considerations. The homogeneous and isotropic state (0| 
admits gapless excitations tt, 

tp — flt + TT, (5) 

the superfluid phonons. At finite temperature, these will 
get excited, will reach thermodynamic equilibrium, and 
will form a thermal bath. If one now considers out-of- 
cqulibrium perturbations, with frequencies lower than 
the phonon inverse mean free time, and with wavelengths 
longer than their mean free path, this phonon bath will be 
describable by ordinary hydrodynamics. One thus has a 
normal fluid made up of phonons, moving in the presence 
of, and interacting with a background superfluid — which 
can also "move" of course. (As we stressed already, from 
the field theory viewpoint the hydrodynamical interpre- 
tation of the superfluid dynamics is somewhat acciden- 
tal. By a moving superfluid we mean a field configuration 
different from (0J.) On the other hand, at much higher 
temperatures, above a critical temperature T c , the spon- 
taneous breaking will be gone, superfluidity will be lost, 
and the dynamics will be those of an ordinary fluid — 
like for any other system. It is thus natural to postulate 
that at any temperature between zero and T c , the system 
be made up of two fluids — a 'super' one and a 'normal' 
one, and that the density of the former decreases mono- 
tonically with temperature, dropping to zero at T c and 
staying zero thereafter. This is, in essence, Landau's two- 
fluid model. 

Taking for granted these physical facts, we seek a low- 
energy field theoretical description of such a two-fluid 
system. We should first isolate the low-energy degrees of 
freedom and the symmetries acting on them. We should 
then construct the most generic Lagrangian involving the 
former and compatible with the latter, organized as a 
derivative expansion. The degree of freedom parameter- 
izing the superfluid component is the scalar ip(x,i) we 
have been considering so far. As we argued, it enjoys 
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a shift simmetry — eq. (JXJ) . The normal component has 
the usual long- wavelength degrees of freedom of ordinary 
hydrodynamics, which for field theoretical purposes are 
conveniently parameterized by three scalar fields cf) 1 (x, t) , 
with I = 1, 2, 3 d These should be thought of as the 
comoving (or 'Lagrangian') coordinates of the fluid el- 
ement occupying physical (or 'Eulerian') position x at 
time t. Their dynamics should be invariant under the 
following internal symmetries 0, Q : 

4>' ->• cf) 1 + a 1 , a 1 = const (6) 

(P 1 -> R'j <f> J , i? G 50(3) (7) 

^-►W). det|^ = l. (8) 

The first two symmetries correspond to the physical 
homogeneity and isotropy of the fluid's internal space. 
The last symmetry is what distinguishes a fluid from 
an isotropic solid — the dynamics' insensitivity to trans- 
verse deformations of the system. The extra symmetry 
we should also impose is Poincare invariance — which is 
straightforward to implement given that we are dealing 
with scalar fields. 

Given the shift-symmetries ([1]) and ([5]), our fields 
should enter the Lagrangian with at least one derivative 
acting on each of them. Therefore, at lowest order in 
the derivative expansion, the Lagrangian is made up of 
d^tp an( i df^cj) 1 . Moreover, in order to obey the volume- 
preserving diff symmetry (JSJ), the fluid variables </> 7 have 
to appear in the combination 

With the two vectors d^ip and J M , we can construct 
three scalar quantities — the three independent scalar 
products — compatible with all the symmetries. In fact, 
it turns out to be slightly more convenient to first factor 
out J^'s norm: 

b = yf-J^J* = ^det d^d^ J (10) 

(J p is time-like — hence the minus sign.) The resulting 
normalized four- vector, 

«* = \j» , (11) 

is nothing but the fluid's four-velocity: it is normalized 
to —1, and, given eq. (9), it is a vector field along which 
comoving coordinates do not change: 

u^ z = 0, 7=1,2,3. (12) 

These properties define the fluid velocity field. The three 
invariant scalars we can construct at this order in deriva- 
tives therefore are 

b, X = d^d^, y = u»d^. (13) 



The low-energy effective Lagrangian should be a generic 
function of them: 

C = F(b,X,y). (14) 

Our claim is that this field theory encodes the infrared 
dynamics of finite-temperature relativistic superfluids. 

Before checking this claim against more standard for- 
mulations of superfluid dynamics, we want to stress the 
economy of our approach, and the conceptual advantages 
that such economy offers. The degrees of freedom are 
parameterized by four scalar fields, with several internal 
symmetries. At low-energies, such symmetries are pow- 
erful enough to isolate three scalar combinations of our 
fields as the only possible invariants. The Lagrangian is 
a function of these three quantities. Everything else fol- 
lows straightforwardly from the Lagrangian: the equa- 
tions of motion (upon varying with respect to the fields) , 
the stress-energy tensor (upon varying with respect to 
the metric), the conserved U(l) current (via Noether's 
theorem), the spectrum of excitations (upon expanding 
at second order in small perturbations), the interactions 
of such excitations (via standard perturbation theory). 

Moreover, we want to stress that unlike for previ- 
ous Lagrangian approaches to fluid-dynamics, our La- 
grangian is really a Lagrangian in the most ordinary 
sense. The equations of motion are obtained via the stan- 
dard variational principle — we vary the fields while keep- 
ing their values at infinity fixed, without further restric- 
tions on the variations. All fields have non- vanishing con- 
jugate momenta, and, consequently, none of them plays 
the role of a Lagrange multiplier. This is achieved by 
treating the positions of the fluid elements as our dy- 
namical variables, rather than bypassing them and try- 
ing to construct an action directly for apparently more 
natural degrees of freedom like the velocity and the den- 
sity fields. In our approach these are derived quantities, 
precisely like in Lagrangian mechanics the velocity q(t) 
is a derived quantity, while it is the coordinate q(t) that 
plays the role of the fundamental 'field' entering the ac- 
tion principle. 

III. FIRST AND SECOND SOUND 

As we will see, because of the unconventional language 
that we are using, the most laborious part of our anal- 
ysis will be matching our results to the standard ones. 
However, to convince the reader of the usefulness of our 
approach, we can show straightforwardly that the spec- 
trum of excitations is composed of two propagating lon- 
gitudinal sound modes, as befits a superfluid at finite 
temperature. 

Consider an isotropic and homogeneous background 
configuration, where ip takes the form (QJ, and the nor- 
mal fluid part is at rest and homogeneous in such a ref- 
erence frame. This means that the comoving coordinates 
(p 1 are aligned with the physical ones, (j) 1 oc x 1 . The pro- 
portionality constant measures the compression level of 
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such an homogeneous state, and depends, for instance, 
on the applied external pressure. Now consider small 
perturbations about such a background: 

1I>(X) =VO-(t+ 7T°(X)) , <?{x) = bl /3 ■ (X 1 + TT^X)) . 

The proportionality constants have been chosen for no- 
tational convenience, since 

y = y + 0(tt) , b = b + 0(tt) . (16) 

To study the free propagation of small perturbations, one 
has to expand the Lagrangian at quadratic order in tt° 
and 7T. After a straightforward Taylor expansion, we get 

C~\[K n tt 2 -GW(V -if) 2 ] 

+ M (V ■ 7r)7T° (17) 
where the various coefficients are defined as 



(F v yo - F b b ) 



-Fbbb 2 , 



K N 
Gn 

K s = (F yy - 2F x )y 2 
G s - ~2F x y 2 

M = F b yb ya - F y y - 2F bx b yl 



AF yX yl + AF xx y 4 



(18) 
(19) 
(20) 
(21) 
(22) 



On the right-hand-sides, the subscript b, y, and X stand 
for partial derivatives of F — all computed at the back- 
ground values b — b , y — y 0l and X = —y 2 ,. On the 
left-hand-sides, K, G, N, S, and M stand for 'kinetic', 
'gradient', 'normal', 'super', and 'mixing', respectively. 

Despite the apparent complexity of these coefficients, 
the structure of the quadratic Lagrangian (|17|) is ex- 
tremely simple. We can see right away that, precisely 
like for an ordinary fluid, the transverse (i.e., divergence- 
less) component of n does not propagate. It enters the 
quadratic Lagrangian as 



(23) 



that is, there is no gradient energy for it, nor does it get 
one from mixing with ir° — this is a direct consequence 
of the volume-preserving diff symmetry jSJ. As a result, 
7Tt does not feature propagating wave solutions. (For a 
more in-depth discussion about this point, we refer the 
reader to Q.) 

On the other hand, 7?'s longitudinal component and 
7T° have standard kinetic and gradient energies. More- 
over, there is a term mixing them — that weighed by M. 
Given that F is a generic function — we will see that it 
is related to the equation of state, which can be fairly 
generic — we should think of all the coefficients (IT%l) - (f2"2")) 
as independent, since they involve different combinations 
of derivatives of F. As a result, barring accidental can- 
cellations, we expect that the diagonalization of the ttl- 
7T° system yield two non-degenerate propagating modes, 



with eigenfrequencies determined by the secular equation 

detf X <"^ 2 „ Y Uk rh *)=0- (24) 
V Mujk K s u l - G s k z J v ; 

We thus have two independent kinds of sound waves, 
propagating non-dispersively (i.e., with uj oc k) at dif- 
ferent speeds, and corresponding to different linear com- 
binations of a longitudinal deformation ttl of the nor- 
mal fluid component, and of a perturbation ir° of the 
supcrfluid one. These are precisely the qualitative prop- 
erties of first and second sound. We will see that our 
description works quantitively as well — for instance the 
propagation speeds match precisely those computed via 
Landau's hydrodynamic theory — but before being able 
to do so, we have to find the dictionary connecting our 
field-theoretical language to hydrodynamic and thermo- 
dynamic variables. 



IV. HYDRODYNAMICS AND 
THERMODYNAMICS 

We will follow the same logic as in [10(. We start by 
computing the stress-energy tensor associated with the 
Lagrangian (|14[) . Let's consider an infinitesimal variation 
of the spacetime metric, g^ v — rf v + 5g^ v . We should 
compute the associated variation of the action, when we 
keep all our field variables ip, (ft 1 constant. The metric 
enters the action in several ways. First of all, there is the 
overall \J—g upfront. Second, 6, X, and y all depend on 
the metric. For b, from the r.h.s of eq. (|10|) we have 



5b = \bBj} d^dvtf Sg^ 



where Bjj is the inverse of 



B IJ = d^d^ , 
and we made use of the identity 

b 2 = det B IJ . 
For X, we have simply 

5X = d^d^Sg^ 



(25) 



(26) 



(27) 



(28) 



For y, the metric enters through u M = J^/b, which de- 
pends on it both because of 6, and because of the e- 
symbol in J M , e M "^ 7 oc Collecting all these contri- 

butions we get 



ss 



flv ~ Sg^ 

= (F y y - F h b)Bj}d tl tfd v <t> J + (F — F y y) n ^ 
- 2Fxd^d u i> . 

Using the identity 

Bjjd^cj) 1 d„(f> J = rj^u + u^U v , 



(29) 



(30) 
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we can rewrite this as 

T^ v = (F y y - F b b)u^u v + (F — F b b) V ^ - 2F x d^d v ^ . 

(31) 

Next, we compute the current associated with the 
spontaneously broken U(l) symmetry ([1]). The La- 
grangian depends on ip only through d^tp. In such a 
simple case, the shift Noether current reduces to 

= dC 

= F y u" + 2F X d^i) . (32) 

Now, to match our field-theory quantities to the more 
standard hydrodynamical and thermodynamical ones, we 
have to extract the energy density p, the pressure p, and 
the particle number density n, from T M „ and above, 
and then to impose the usual thermodynamics identi- 
ties. To extract p and n we have to pick a local reference 
frame. Since we have two fluids in relative motion, there 
are two natural reference frames, which means that nei- 
ther is really the preferred one. Ultimately it does not 
matter which frame we pick as long as we apply ther- 
modynamics consistently in that frame. We choose to 
move with the normal fluid component, along u^, also 
because — as already emphasized — the interpretation of 
d^ip as a superfluid velocity field is perhaps possible, but 
not necessary. We have 

p = T» v u„u v = F y y-F- 2F x y 2 (33) 
n = -f\ t =F y - 2F x y (34) 

(the minus sign in the definition of n takes care of the 
negative normalization of u^.) As to the pressure, unlike 
for a single fluid, there is in general no frame in which 
the stress-tensor is isotropic. This makes the definition 
of pressure less obvious than usually. However, is 
still the sum of an isotropic tensor (oc r]^ v ) and of rank- 
one tensors, one proportional to u^ L u v , and another one 
proportional to d^tpd v tj). This decomposition is unam- 
biguous, and invites identifying p with the coefficient of 
r\y, v . This means that in any reference frame, and for non 
relativistic fluid motions, the stress tensor Tjj will be the 
sum of pSij and of two anisotropic terms, proportional 
to v l N v° N and to where v x and vs are the normal 

fluid and superfluid velocity fields. This is the standard 
definition of pressure for non-relativistic superfluids, and 
from the thermodynamical viewpoint it corresponds to 
the derivative of the total energy w.r.t. the volume, taken 
at constant total entropy, charge, and relative momentum 
0. We are thus led to set 

p = F-F b b. (35) 

We now impose the thermodynamic identities. We 
start with 

p + p = Ts + pn . (36) 



Plugging in the expressions for p, p, and n derived above 
we get 

(F y - 2F x y)y - F b b = Ts + p(F y - 2F x y) (37) 

This suggests the identifications 

p = y, s = 6, T = -F b , (38) 

which are consistent with Landau's idea that the entropy 
is carried by the normal fluid component only: b does not 
depend on the superfluid field variable ip, but only on the 
normal fluid ones, (p 1 . Moreover, as emphasized in [To| . 
with these definitions the entropy current coincides with 
our J M = 6u M , which is identically conserved by virtue 
of its definition (J9j) . There is course an ambiguity in the 
overall normalization of s and T — we can multiply s by a 
constant and divide T by the same constant without af- 
fecting the thermodynamic identities. This just amounts 
to changing ks, that is, to a change of units for temper- 
ature. 

Then, we should impose one of the differential thermo- 
dynamic identities — for instance 

dp = s dT + n dp + . . . (39) 

The dots stand for an additional term that takes into 
account that here the number of independent thermo- 
dynamic variables is three, rather than the usual two. 
This is evident from our ficld-thcoretical description: the 
dynamics encoded by the Lagrangian (|14[) involve three 
independent scalar quantities. We have already identi- 
fied b with the entropy density, and y with the chemical 
potential. We are still missing the thermodynamic inter- 
pretation of X. In more standard treatments of (non- 
relativistic) superfluid dynamics, the extra independent 
variable is the relative velocity — or the relative momen- 
tum density — between the superfluid and normal fluid 
parts. We will recover this interpretation below. 

If we compute the differential of (1331) in terms of dT = 
—dF b , dp = dy, and dX, we get 

dp = -b dF b + F y dy + F x dX 

= sdT + ndp + F x {dX + 2ydy) . (40) 

The combination in parentheses is most easily written in 
terms of the components of d^tp orthogonal to u^: 

drf = -u M v + b, e = W + ^u»)d^j. (41) 

being orthogonal to u'*, can be thought of as a three- 
vector £. We have X = —y 2 + £ 2 , so that 

dp = sdT + n dp + 2F X £ d£ , (42) 

with £ = 

At the level of hydrodynamics and thermodynamics, 
our description of the system is manifestly equivalent to 
that of Son et al. [TH, H3| , which in turn is equivalent to 
those of Israel Q , of Khalatnikov and Lebedev Q , and 
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of Carter Q (see [1, 0] for the equivalence of these differ- 
ent classic approaches). To see this, apart from obvious 
changes in the notation, one only needs to take into ac- 
count that n and p here stand for the total number and 
energy densities. On the other hand, in [l3|, [l4| n stands 
for the normal component's number density — our F y in 
eq. (I34[) — while e is defined via e + p = Ts + \in (20| . 

To summarize, we see from eq. (|42p that pressure is 
naturally a function of T , [i, and £. Such a function 
defines the equation of state of our superfluid. Once it is 
given, one can construct the corresponding Lagrangian 
(fill from eq. l[35]). that is 



C = F = p + F b b 

= P ~ T df> 
and @2J}. 



where we used eqs. 

then be expressed in terms of b = s = y 



(43) 
(44) 

The result should 
/i, and 



X = ? 



We will see an explicit example of this 



procedure below. 



Relative motion 



One of the surprising features of superfluids at finite 
temperature is the possibility of having, in thermal equi- 
librium, a relative motion between the two components. 
The two fluids can freely flow through each, provided 
their motion is homogeneous and happens at constant 
speed. In our language, a necessary condition for this 
to be possible is that our field equations admit station- 
ary and homogeneous solutions where the reference frame 
defined by the superfluid — which is associated with the 
vector field d^ip — is at not rest with respect to the nor- 
mal fluid's one — which is associated with the vector that 
we call m m . In other words, there should be solutions of 
the form, 

7 (x, t) — box 1 , ip(x, t) = j/o t + £o • x , (45) 

where bo, yo, and £o are constant. In particular, £o is 
the same variable we introduced in the last section, and 
it quantifies the misalignment between the two reference 
frames: 



(1,0), d^ = (y ,£ ) 



(46) 



It is easy to see that all configurations of the form 
(1431) solve our field equations. The reason is that the 
Lagrangian (|I4[) only depends on first derivatives of the 
fields. As a result, the field equations all take the 
schematic form 



dG(d<p I ,dip) = 



(47) 



Such equations are obviously obeyed by the field config- 
uration (1451). which has constant first derivatives. 



The existence of such solutions is only a necessary con- 
dition for relative motion's being compatible with ther- 
mal equilibrium. This is because our field theory neglects 
dissipative effects, which are clearly crucial for reaching 
equlibrium. What our solutions show is the absence of 
conservative forces between the two fluids in relative mo- 
tion. However as we argued above, for our field theory 
to describe a thermal system, this has to be character- 
ized by three thermodynamic variables, one of which is 
precisely £. In other words, the equation of state — which 
is a statement about the system at equilibrium — depends 
on the state of relative motion. This means that relative 
motion is compatible with thermal equilibrium. 



V. THE HIGH-TEMPERATURE LIMIT 

The high temperature behavior of a superfluid is char- 
acterized by the disappearance of superfluidity at some 
critical temperature T c . Right beneath T c , the superfluid 
component is extremely small, and we can work pertur- 
batively in it. That the superfluid component is small, 
concretely means that its contributions to the stress- 
energy tensor and to the charge current are small. From 
eqs. (I3T1) . (|3^|) we see that this corresponds to having 
small Fx' in the limit in which Fx vanishes, the stress- 
energy tensor and the current reduce to those of an or- 
dinary fluid flowing along , which for us is the normal 
component's four-velocity field. We are thus led to the 
conclusion that the dynamics of a superfluid close to the 
transition temperature are captured by the behavior of 
our field theory in a region of field space where Fx is 
small. In fact, at T = T c and above, the superfluid com- 
ponent disappears completely, and one is left with an 
ordinary fluid carrying a conserved charge. In our La- 
grangian formalism, this is described by [10j 



•<T>T C 



F(b,y) 



(48) 



This means that above T c there is no dependence on X 
whatsoever. As pointed out in fiof . this is equivalent to 
imposing the generalized shift-symmetry 



(49) 



for generic /, under which b and y are invariant, but X is 
not. Therefore, even though it looks "unnatural" to have 
a whole region of (b, y, X)-space where F does not depend 
on X, such a region is in fact one of enhanced symmetry. 
As a result, its existence is protected by symmetry. 

Let's consider the case described by (f4"5)l first. For 
small perturbations about the equilibrium configuration, 
the quadratic Lagrangian (|I7|) simplifies to 

C ~ \ [(F y y - F b b )ir 2 + F bb b 2 (V • t?) 2 ] 

+ h F vvVl (^°) 2 + (Fbyboyv - F v y )(V ■ tt)tt° . (50) 

Notice that now there is no gradient energy for tt°, but 
there is still a mixing between ir° and ttl- The La- 
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grangian is easily diagonalized via the field redefinition 
*° = ^-%S^--i), (51) 

7?L = TTl , (52) 



where 7?i is purely longitudinal. Neglecting 7tV — which 
does not propagate — we get 



,lp 2-2 



F yv F bb b z a -(F by b -F y )- _ X21 

(53) 



We see that 7r 2 — like ttt — does not have a gradient en- 
ergy term. As a consequence, it does not feature wave 
solutions. On the other hand, 7?i has wave solutions prop- 
agating at (squared) speed 



(F by b - F y f - F vy F bb bl 
Fyy(F y yo - F b b ) 



(54) 



Equivalently, when all derivatives of F with respect to X 
vanish identically, the secular equation (124[) has a non- 
trivial solution, ujf = cfk 2 , corresponding to waves of ffi, 
and a trivial one, lo 2 = 0, corresponding to 7r2's trivial 
dynamics. 

We thus see that above T c there is only one propagat- 
ing sound mode, as befits an ordinary fluid. It is some- 
what non-trivial to check that the propagation speed we 
found — eq. (1541) — matches the standard hydrodynamic 
expression 



dp 
dp 



(55) 



S.N 



where S and N are the total entropy and charge. The 
reason is that our Lagrangian is naturally a function of 
b = s and y = p, rather than of s and n — and so are 
p and p in our formalism. A similar complication was 
encountered in and it was overcome by using ther- 
modynamic identities. Here, we can use the following 
trick. In (|55|) , the constraint of being at constant S and 
at constant N gives a one-to-one relation between dp and 
dp. In our field theory, the same relation can be obtained 
via the equations of motion: the dynamics of our system 
conserve the total entropy by construction (since they fol- 
low from a conservative Lagrangian) , and the total charge 
because of a U(l) symmetry. Whatever relation we get 
between dp and dp by enforcing the equations of motion, 
is thus going to be equivalent to that enforced in (|55|) 
plj . Charge current conservation reads 



WO 



o . 



(56) 



(we are setting Fx = 0). By using it M = bJ**, and the 
conservation identity for J M (which is equivalent to the 
conservation of the entropy current), we get 

J^(F y /b)=0 (57) 

That is, (Fy/b) is conserved along the flow, or, 



bdF, 



F y db 



(along the flow) 



(58) 



We thus see that the equations of motion imply a rela- 
tion between dF y and db. This in turn translates into a 
relation between dp and dp. To see this, let's first express 
dF y in terms of dy and db: 

dFy = Fyydy + Fy b db . (59) 

This implies a relation between dy and db along the flow; 

bFyy dy = (F y - F yb b)db (along the flow) . (60) 
We thus get 



dp 
dp 



S,N 



dp 
dp 

d(F - F b b) 



flow 



d(F v y 
(Fbyb- 



-F) 

Fy) 



flow 



FyyFbbb 



Fyy{F v y-F b b) 



(61) 



which is precisely eq. (|54"1) . 

At temperatures slightly below T c , F acquires a weak 
X-dependence. Derivatives of F with respect to X are 
now non-vanishing, but can be treated as small. In par- 
ticular, in the coefficients entering the quadratic action 
for small excitations, eqs. ([T% |) - (|!21?]) . we can ignore these 
X-derivatives everywhere except for Gs- The reason is 
that Gs vanishes when these derivatives do, and as a con- 
sequence the secular equation becomes degenerate — 
i.e., one of the solutions becomes trivial, and one mode 
stops propagating. When Gs ~ Fx is small but non- 
vanishing, this degenerate mode will start propagating, 
with a frequency w| = 0(Fx)k 2 . Keeping into account 
the X-derivatives of F for the other coefficients would 
simply shift the two eigenfrequencies by relatively negli- 
gible amounts. That is, schematically: 



,? = ( C ? + 0(F x ))fc 2 

4 = 0(F x )(l + 0(F x ))k 2 



(62) 
(63) 



To compute the speed of the new propagating mode at 
leading order, we just have to keep Fx in Gs and drop all 
other X-derivatives from the other coefficients. The rele- 
vant quadratic Lagrangian is then eq. (|50|) supplemented 
with a small gradient energy for 7r°: 

C ~ | [(F y y - F b b )ir 2 + F bb b 2 (V • t?) 2 ] 

+ ±[Fy y y 2 (n°) 2 + 2F x y 2 (^ ) 2 ] 

+ (Fbyboya - F y y )(V ■ ^)tV° . (64) 

The near vanishing eigenfrequency of the secular equa- 
tion becomes 

2 _ G N Gs , 2 



2 M 2 + G N K S 

2F bb b 2 y 2 



(F by b Q y - F y y ) 2 - F bb F yy b 2 y 2 



F x k 2 



(65) 
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It is straightforward but somewhat tedious to check that 
in the non-relativistic limit, this reduces precisely to the 
second sound dispersion law as predicted in Landau's hy- 
drodynamic theory for superfiuids near 

For the same reasons as above, for finite but small Fx 
the eigenmodes are still 7?i and tt2 as defined in ([5Tjl . plus 
small 0(Fx) corrections which we can neglect at leading 
order. 



VI. THE LOW-TEMPERATURE LIMIT 

Things are much neater at very low temperatures. At 
exactly zero-temperature, we just have the supcrfiuid 
component, and the Lagrangian depends on X only, as 
in eq. ©: 

C T =o = F(X) . (66) 
All kinetic coefficients (IT8"1) - ([2"2"|) vanish except for 

K S = ^Fxxvt - IFxyl , G s = -2F x y 2 . (67) 
The propagation speed of ir° waves is thus 
2 F x - 2F XX yi _ dp 



F 



x 



dp ' 



(68) 



where the subscript 'I' stands for 'first sound'. Notice 
that dp/ dp here is unambiguous, since both p and p de- 
pend on X only. 

At very small but non-vanishing temperatures, these 
phonons get excited a make up a normal fluid. This 
changes the equation of state, and consequently the La- 
grangian gets slightly modified 



£ = F(X) + 5F(b,y,X) 



(69) 



The entropy density is dominated by the phonon contri- 
bution, which scales as s ~ T 3 . Since for us s = b and 
T = —5F b , we expect that at low temperatures 



SF = b 4 ^f( yi X) b^O. 



(70) 



We will see below that this is indeed the case. All kinetic 
coefficients pl)-(|22]l get modified by C(6 4 / 3 ) = 0(T 4 ) 
corrections. However, for the same reasons as those an- 
alyzed above, we need to keep track of these corrections 
only for those coefficients that would vanish otherwise, 
which are Kjy, Gn, M. Moreover, since M contributes 
an off-diagonal term to the kinetic matrix, it appears 
quadratically in the secular equation, (|24p. while Kn and 
Gn appear linearly. M can thus be neglected at leading 
order. The low-temperature quadratic Lagrangian then 



C ~ \ [(SF y yo - 5F b b )T? 2 + SF bb b 2 Q (V ■ if) 2 ] 
+ \ [^Fxxy\ ~ 2Fxyl){^f + 2Fxy 2 (Vn ) 2 



It is already in a diagonal form — the two sound modes 
coincide with ir° and ttl , and have propagating speeds c\ 
as given in f[68]) and 
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SF bb bl 



SF y yo - 5F b b 



(71) 



respectively. 

Notice that in this limit the second sound corresponds 
to longitudinal perturbations of the phonon gas. Even 
though such a gas fades away as T approaches zero, the 
second sound speed approaches a finite value. Indeed, 
for non-relativistic superfiuids Landau predicted that it 
should approach the limiting value c| — > cf/3. This re- 
sult has been generalized to relativistic superfiuids in 
[l6|, where the same limiting value was found. To see 
whether our expression for c 2 , reproduces this result, we 
need the low-temperature equation of state. Fortunately 
the finite-temperature part of such an equation of state — 
which is associated with our SF — is model- independent. 
The reason is that it is determined purely by the low- 
temperature behavior of the phonon gas. Since phonons 
are derivatively coupled to each-other, at low tempera- 
tures they behave as free particles propagating at speed 
ci . c\ is thus the only parameter of the underlying zero- 
temperature superfluid entering the leading order finite 
temperature correction to the equation of state. 

The low-temperature equation of state for relativistic 
superfiuids has been computed in fig . In our notation, 
it relates the pressure to b, y, and X: 



p~F(X) 
Using (|35|) we get 



SF(b,y,X) 



-0 

ci V 



1/3 



^1 
Cl 



(l + (!-"*>£)' 



1/3 



(72) 



(73) 



which, indeed, is of the form (|70[) . Evaluating the 
derivatives entering (|TTj) at the background values bo, yo, 
Xo = -Vo, we get 



3 C 1 7 



(74) 



as expected. 



VII. SOUND WAVE SCATTERING OFF A 
SUPERFLUID VORTEX 

So far we have been preoccupied with trying to dig- 
nify our field-theoretical approach by making it repro- 
duce well known results about superfiuids. Here instead, 
we make it compute something new. One of the main ad- 
vantages of a field-theoretical Lagrangian description, is 
that it offers a systematic and straightforward approach 
to perturbation theory, e.g. for the study of scattering 
processes. With our superfluid theory we can consider 
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the scattering of sound waves colliding with each other, 
along the lines of @. Perhaps more interesting is the 
scattering process involving a sound wave and a super- 
fluid vortex, which has been claimed to be a powerful 
probe of 'quantum vorticity' [17] . A similar process has 
been considered in [12| for supersolids. 

Consider a superfiuid vortex configuration in an oth- 
erwise homogeneous, still superfiuid: 

tp(x,t) = y t + ip v , ip v (x,t) =Cip (75) 
(/) I (x,t) = bl /3 x I . (76) 

We are assuming that the vortex is aligned with the z- 
axis, and we are denoting the azimuthal angle by if. C is 
a constant that determines the periodicity of ip: since ip 
is the Goldstone boson for a spontaneously broken U(l) 
symmetry, it should be periodic, ip ~ ip + 2ttC. Without 
loss of generality, we can set 



C = 1 



(77) 



This just fixes the normalization of ip, and of the associ- 
ated charge. 

Before considering the scattering of sound waves off 
the vortex, we have to make sure that the vortex con- 
figuration is a solution. Since the presence of the vortex 
perturbs the spatial derivatives of ip only, y and b are 
unaffected by it: 



y = u^dutp = 2/0 



b = b 



On the other hand, X is perturbed substantially: 
X = (di'f = -y 2 + \ . 



(78) 



(79) 



This suggests that if the Lagrangian did not depend on 
X, then ([751 ITB")) would be a solution. To check that this 
is indeed correct, consider the Lagrangian (|48|) . As we 
mentioned, it enjoys a generalized shift-symmetry, ip — ¥ 
ip + /((f) 1 )- Now, the vortex configuration ([751 1761 can 
be obtained by applying one such transformation to the 
background configuration 



with 



ip {x, t) = y t , (f>l (x, t) = bl^x 1 



= arctan^o/^ = (p. 



(80) 



(81) 



Therefore, since the background (1801) solves the equations 
of motion, so does the vortex configuration. 

Now, the Lagrangian ([48]) is the correct description of 
the system at or above the critical temperature. At tem- 
peratures right below it, the Jf-dependence of F can still 
be treated as small, which means that the true vortex 
solution will be slightly different from ([751 ITS"]) . At lead- 
ing order we can neglect this small difference and stick 
with ([751 176[) . At significantly lower temperatures, there 



will still be vortex solutions, but they will differ signifi- 
cantly from ([751 175)) . In particular, we expect the nor- 
mal component — parametrized by the 7 's — to be signif- 
icantly perturbed. The study of these more complicated 
vortices is beyond the scope of this paper. We focus on 
temperatures right beneath the critical one, where (J75J 
[76]) is the leading order vortex solution. 

To study the scattering of a sound wave off the vortex, 
we need the quadratic Lagrangian for small perturbations 
in the presence of the vortex: 

^ = yo(< + 7r°) + Vv, ¥ = bl /3 ■ (x 1 + tt 1 ) . (82) 

At leading order in the X-derivatives of F, we get ([M]) for 
the vortex-independent part of the quadratic Lagrangian, 
and, after some integration by parts, 

dot ~ [(F y - F yb b )(V • tt) - F yy y tt°]tt • Vipv 

- \F y riV (OA - d k di)ip v (83) 

for the vortex-sound wave interaction part. The Taylor 
expansion leading to this expression is straightforward. 
We have omitted O(diipv) 2 terms, since their relative 
importance in a scattering process is suppressed at low 
enough sound-wave frequencies (lower than the chemical 
potential). The anti-symmetrized second derivative ap- 
pearing in the second line does not vanish, due to the 
non-trivial topology of the vortex solution: 



{d x d y - d y d x )ip v = (2tt) 5 2 {x, y) 



(84) 



Precisely the same 'topological interaction' term was 
found in [13] in the context of supersolids. 

The interaction Lagrangian further simplifies if 
we express it directly in terms of the propagating 
eigenmodes — first and second sound. At leading order 
in Fx they are given in (|5ip. We get simply: 



(85) 



Ant — ~F yy yo K 2 7Tl ■ Vtjjy 

-tt F y (tti x 7?i) • zS 2 (x, y) 



The first line describes the conversion of second sound 
into first sound, and vice versa. The second line describes 
the elastic scattering of first sound off the vortex. At 
this order in Fx, there is no elastic second sound- vortex 
scattering — but this may be deceiving: the second sound 
speed also vanishes in the Fx — ► limit. The incoming 
and outgoing states' propagation speeds enter the scat- 
tering cross section in several ways, e.g. via the phase- 
space element. As a consequence, it is possible that the 
second sound-second sound scattering cross section ap- 
proach a finite limit for Fx 



17] claims that 



0. Indeed, 

elastic second sound scattering by a superfiuid vortex is 
the dominant process for superfiuid helium. Analyzing 
this possibility is beyond the scope of this paper, and we 
leave it for future work. 

For simplicity, let's restrict to scattering processes 
where the incoming wave travels orthogonally to the 
vortex, k ■ z — 0. Since the vortex does not break 
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translations along z, the z-component of the momen- 
tum is conserved. This implies that the outgoing wave 
is also going to propagate perpendicularly to the vor- 
tex, k! ■ z — 0. This effectively reduces the problem to 
a two-dimensional one, taking place in the x-y plane. 
Moreover, energy is also conserved, w' = w, because the 
vortex does not break time translations. To compute the 
associated scattering amplitudes, we apply the standard 
Feynman rules associated with the so-called relativistic 
normalization for the single-particle states. As reviewed 
in Q , this approach is consistent — and convenient — even 
when Lorentz-invariance is not there (in our case, it is 
spontaneously broken by the background). We have to 
take into account the non-canonical normalization of our 
first and second sound fields — see eq. (|64|) . To simplify 
the notation, let's denote the combination (F y yo — Ff,bo) 
by w. The relevant (two-dimensional) S'-matrix elements 
then are 



iMi^i = -2ni ^- -ujsmt 
i— ?i w 

iM^ 2 = iM 2 _,! 



(86) 



c 2 uj sm9- (2ir)5{u' - u) , (87) 

where 8 is the scattering angle. Since the vortex breaks 
translations in the x-y plane, there is no (^-function en- 
forcing two-dimensional momentum conservation. For 
the 1 — > 2 and the 2 — > 1 conversion processes, we have 
kept the leading order in c 2 /ci <C 1. 

The associated cross-sections per unit vortex length 
are 

dlJa-fb 111.. 1 2 1 



d6 dz 2lu c a 



M a - 



Attc 2 



where c a ^ are the initial and final states' propagation 
speeds, Ai a ^b are the S'-matrix elements M a ^,b without 
the (2tt) S(. . . ) factors, and the last factor comes from 
the final state's phase space element: 

cP-V 1 

(2ir)5(u/-u) - 



dll h = 



d9 1 



We get 



(2tt) 2 2u> 



da i. 



IT- (89) 



4tt c\ 



7T 



d9 dz 
d6 dz 

d(72^1 

d6 dz 



F 2 

v 

2 w 2 



1 



it F 
2 w c\ 
C2 dai^ 2 
ci 



dd dz 



(90) 
(91) 
(92) 



We can express these formulae in terms of thermodynam- 
ical quantities. To this end, recall that for near vanishing 
F x , we have (see sect. IIVI) 



Fy — 1 
We thus get 



w 



F y y-F b b~p + p, F y 



dn 
dfi, 



d<7i 



d9 dz 

d<Ji^ 2 
d0 dz 

dP2^\ 

d6 dz 



2(p + p) 2 
7T (dn/dfi) t 
2 

Cl 



(p+p) 

do\^ 2 
dd dz 



Cl 



(93) 

(94) 
(95) 
(96) 



These results are fully relativistic. In the non-relativistic 
limit, (Ti_).i reduces precisely to the elastic phonon- vortex 
cross-section computed by Son for supersolids 

A final comment is in order: to compute these scat- 
tering cross sections we have found it convenient to use 
a quantum mechanical language. However our results 
are entirely classical, and hold for classical sound waves 
without any modifications. As usual for the scattering of 
waves off a 'potential' or an external source — the vortex, 
in our case — one can use tree-level quantum perturbation 
theory to compute the classical cross-section, defined as 
the ratio between the scattered power and the incoming 
energy flux. 



VIII. DISCUSSION 



We have presented a low-energy, long-distance effec- 
tive field theory for finite-temperature relativistic super- 
fluids. It is equivalent to more standard hydrodynami- 
cal approaches, yet it offers — we believe — a conceptually 
clearer and technically simpler framework. As empha- 
sized in [10] for ordinary fluids, it may offer a particu- 
larly convenient setup to organize the derivative expan- 
sion. Besides, certain properties may be invisible — or 
at least very well hidden — at the level of hydrodynam- 
ical equations, yet manifest at the level of the action, 
like for instance the unitarity constraints of [l8j]. More- 
over, the Lagrangian formalism is invaluable for pertur- 
bation theory, e.g. in the systematic study of scattering 
processes. As an example we have computed the cross- 
section for sound wave-vortex scattering. The computa- 
tion is straightforward, and the result easily expressed in 
terms of hydrodynamical quantities and of the equation 
of state. One of the possible applications of our formal- 
ism we envisage for the near future, is reproducing the 
hydrodynamical effects of quantum anomalies [l9] within 
effective field theory. This would require supplementing 
our field theory with Wess-Zumino-like terms, along the 
lines of EI. 
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